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We examine linear dust acoustic waves (DAWs) in a dusty plasma with strongly correlated dust 
grains, and discuss possibility of a twisted DA vortex beam carrying orbital angular momentum 
(OAM). For our purposes, we use the Boltzmann distributed electron and ion density perturbations, 



the dust continuity and generalized viscoelastic dust momentum equations, and Poisson's equation 



to obtain a dispersion relation for the modified DAWs. The effects of the polarization force, strong 
dust couplings, and dust charge fluctuations on the DAW spectrum are examined. Furthermore, 
we demonstrate that the DAW can propagate as a twisted vortex beam carrying OAM. A twisted 
DA vortex structure can trap and transport dust particles in dusty plasmas. 
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I. INTRODUCTION 



Charged dust grains [l| and dusty plasmas [2[ are ubiquteous in cosmic and astrophysical 
environments 3r|5|, such as interstellar media, molecular dusty clouds, star forming dust 
clouds, Eagle nebula, and supernovae remnants, etc. They are also found in our solar 
system, e.g. in planetary rings systems j6-8|, in interplanetary media jo| due to the presence 
of cometary dust particles, on the Martian surface as dust devils, on the surface of Sun 
and on moon, as well as in the Earth's mesosphere 
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12], in space as charged dust debris 
ij produced when satellites are deducted, aad near space propulsion vehicles Q [for 
future spacecrafts that go nearer to the Sun (such as Solar Orbital and Solar Probe Plus) 
and Lunar Atmosphere and Dust Environment Explorer mission to the moon] for exploring 
composition of dust grains and their role in collective dust-plasma interactions) due to rocket 
exhausts. Charged dust particles, which could be of different sizes (ranging from micron-sized 
to nanometer sized), are naturally formed in industrial processing [15| for nanotechnology 
and in magnetic fusion reactors la] . Furthermore, low-temperature dusty plasmas are 
also produced in laboratory devices [l7] for fundamental studies in a new environment 
(e.g. on board International Space Station for examining the behaviour of dusty plasmas 
under microgravity conditions [18|) that does not exist in the usual electron- ion plasma 
without dust. It emerges that dusty plasmas are of broad interdisciplinary interest in physical 
sciences, and share ideas with other fields, e.g. condensed matter physics and astrophysics, 
urthermore, dusty plasmas are also finding applications in medical and biological sciences 



19]. 



A neutral dust particle in plasmas is charged both negative and positive due to a variety 
of physical processes 



201 ] . photo emissions 



20 4251] . including absorption of electrons from the background plasma 
12 1, triboelectric effects 25j, etc. A dusty plasma is usually composed 



of electrons, positive ions, negative or positive dust grains, and neutral atoms. When the 
interaction potential energy (= Z%e 2 /d, where Zd is the dust charge state, e the magnitude 
of the electron charge, and d the inter-grain spacing or the Wigner-Seitz radius) between 
two neighboring dust particles is much larger (smaller) than the dust kinetic energy rC^T^, 
where ks is the Boltzmann constant and Td the dust particle temperature, the dusty plasma 
is in a strongly (weakly) coupled state. 

More than half a century ago, Wuerker et al. 26| showed that electrically charged iron and 
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aluminium particles having diameters of a few micrometers can be contained in a confined 
space by alternating and dc electric fields. Under the three-dimensional focusing alternating 
gradient focusing force and the Coulomb repulsive force, charged iron and aluminium parti- 
cles form a crystalline array, which can be melted and reformed. This seems to be the first 
indication of ordered dust particle structures in an external confined potential. However, 



following Ichimaru's idea 



271 ] of one- component strongly coupled electron system, Ikezi 28 ] 



postulated the solidification of charged dust particles in a dusty plasma |2|, l38M4l| . when 
the dusty plasma T = Zje 2 /(HcbT^ is close to 172. Such values of T can be achieved in low- 
temperature laboratory dusty plasma discharges at room temperatures owing to the large 
a micron-size dust grain would acquire by absorbing electrons from the background plasma. 
The formation of dust Coulomb crystals and ordered dust particle structures have since 
been observed in the sheath region of many laboratory experiments 29|-|34(. The ordered 



dust particle structures are attributed to attractive forces [35|, |36|, |4l| between negative dust 
grains due to ion focusing and ion wakefields in a dusty plasma sheath with streaming ions, 
as well as due to overlapping Debye spheres 37] and dipole-dipole interactions j^j. 

The collective behavior of dusty plasmas involving an ensembles of charged particles was 
recognized through the prediction of the dust acoustic wave (DAW) by Shukla at the 
First Capri Workshop on Dusty Plasmas in May of 1989, where he suggested the existence of 
the nonlinear DAW in the presence of Boltzmann distributed electrons and ions, and massive 
charged dust particles. Shukla's idea was then worked out in the first paper 43] on linear 
and nonlinear DAWs. It must be stressed that there does not exist a counterpart of the 
DAW in an electron-ion plasma without charged dust grains, since the DAW is supported by 
the dust particle inertia, and the restoring force comes from the pressures of the inertialess 
Boltzmann distributed electrons and ions. Thus, similar to the Alfven wave in an electron- 
ion magnetoplasma without dust, the DAW is of fundamental importance in laboratory and 
space plasmas. The DAW is usually excited by an ion streaming instability 45, 1^6 ], and 
has a frequency much smaller than the dust plasma frequency, extending into the infra-sonic 
frequency range when the dust particles are anomalously heated. The low-frequency (of the 
order of 10 Hz) DA fluctuations were first indirectly observed in the experiment of Chu et 



al. 



291 ] prior to dust particle crystallization, and have since been spectacularly observed in 



many laboratory experiments world-wide 



47H49J, and also in the Earth's ionosphere 50]. 



Thus, the existence of the DAW in dusty plasmas have been demonstrated at kinetic levels, 



3 



and the visual images of the dust acoustic wavefronts by naked eyes are possible |47|. 

In this paper, we revisit the linear DAW in a dusty plasma by incorporating the effects 
of dust particle correlations jsi|, the polarization force due to interactions between 
thermal ions and highly charged dust grains, and dust charge fluctuations (DCFs) [53] . A 
linear dispersion relation is derived and analyzed. The underlying physics of the DAW has 
been put on the firm footing. Furthermore, we also discuss the possibility of a twisted dust 
acoustic wave (TDAW) carrying OAM [we also refer it as a dust acoustic vortex (DAV) 
beam]. A TDAW or a DAV beam can, in turn, be used for trapping and transporting 
charged dust grains from one region to another in laboratory and space dusty plasmas. 



II. THEORETICAL CONSIDERATION 



Let us consider an unmagnetized dusty plasma composed of inertialess electrons and ions, 
as well as strongly correlated negative dust particles of uniform sizes. In the presence of 
ultra-low frequency DAWs, with u <C v en , v in <C k 2 V^ e Ti /\u>\, where uj is the wave frequency, 
i/ e „ (ui n ) the electron (ion)-neutral collision frequency, k the wave number, and Vr e (Vn) the 
electron (ion) thermal speed, both electrons and ions obey the Boltzmann law (deduced from 
the balance of the electric force and pressure gradients of the electrons and ions), since they 
can be considered inertialess on the timescale of the DAW period, and rapidly thermalize 
due to collisions with neutrals. Thus, the electron and ion number density perturbations 
(wei.ii < n e0:i0 ) are, respectively, given by [43(| 

n el w n e0 ——, (1) 



and 



ecp 

n a « -n i0 ——, (2) 



where n e0 and n i0 are the unperturbed electron and ion number densities, respectively, <p 
the electrostatic potential of the DAW, and T e (Tj)the electron (ion) temperature. At 
equilibrium, we have the quasi-neutrality condition |43|, |44j, viz. n^o = n e0 + ZdUdo, where 
Z d0 is the average number of electrons residing on a dust grain, and n d0 the unperturbed 
dust number density. 

The dynamics of dust particles in a our dusty plasma is governed by the hydrodynamic 
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equations composed of the dust continuity equation 

dn dx 



at 



+ "doV • v d = 0, 



(3) 



and the generalized dust fluid momentum equation 
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l + Tm dt 



dv d Z do e u d k B T d 

-ttt + Vd^ d (1 - R)\q> H Vn dl 

dt m d p d 

= ^V 2 v d + -^tliv(V-v d ) 

Pd Pd 



(4) 



where n rfl (<^C n d0 ) and v d are the dust number density and dust fluid velocity perturbations, 
respectively, m d the dust mass, p d = n do m d the dust mass density, R = Z d0 e 2 /AksTiXui 
is a parameter determining the effect of the polarization force |52j that arises due to the 
interaction between thermal ions and negative dust grains, p d n do k B T d = P d the effective dust 
thermal pressure, [i d = l + (l/3)u(T) + (T/9)du(T)/dT the compressibility, T = Zj e 2 /dksT d 
the ratio between the dust Coulomb and dust thermal energies, d = (3/4:7m d0 J 1 ' 3 the Wigner- 
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59] 



Seitz radius, u(T) is a measure of the excess internal energy of the system, which reads 
u(T) ~ -( v / 3/2)r 3 / 2 for T < 1 {viz. a liquid-like state), and u(T) = -0.80r + O^r 1 / 4 + 
O.lOr" 1 / 4 — 0.81 in a range 1 < T < 200. The coupling parameter in dusty plasmas including 
the shielding of a negative dust grain by electrons and ions reads T g = Texp(— k), where 
k = a d /\ D , a d is the inter-dust grain spacing, and = \De^Di/(^ 2 D e + ^Di) 1 ^ 2 the effective 
Debye radius of dusty plasmas, with A/? e = (/cBT e /47m e0 e 2 ) 1//2 and Xdi = (ksTi/ 1 A-Knioe 2 ) 1 ! 2 
being the ion and electron Debye radii, respectively. The dust-neutral collision frequency is 



571 ] v dn = (8/3)v^27rm„n n r^VTn/md, where m n is the neutral mass, n n the neutral number 



density, r d the dust grain radius, Vr n = {kBT n /m n ) l l 2 the neutral thermal speed, and T n 
the neutral gas temperature. The visco-elastic properties of the dust fluids are characterized 
by the relaxation time [55|, |56[ r m = [(£ + 4:r]/3)/n d oT d }/ [1 — \i d + 4w(r)/15], involving the 
shear and bulk viscosities r] and £. There are various approaches for calculating rj and £, 
which are widely discussed in the literature 59(. We note that the generalized viscoelastic 



dust momentum Eq. (4) is an extension of Kaw and Sen (51|, by including the effects of the 



polarization force 



521 ] . The viscoelastic momentum equation, similar to Eq. (4), has also 



been used in the study of collective phenomena in fluids 54( and in one-component plasmas 
with strongly correlated electrons 55l. 156]. 

The DA wave potential <fi is obtained from Poisson's equation 



4vre (n el - n iX + Z do n dl + Z dl n d0 ) 



(5) 



where the dust charge perturbation Zdi(<^. Z d o) is determined from the charging equation 

(6) 



d \ „ a d v 2 



Here we have introduced the notations 



and 



V2 = ad 



pe 



UJ 



X 



De 



exp(-?7 e , 



pi 



Y^exp(-r/ e ) + (1+77;)-^- 



(7) 



(8) 



with T) e) i = Z d0 e 2 1 'adhsT^i. Equation (6) reveals that the DCFs cause an adverse phase 
lag between Zdi and <f), which lead to the DAW damping 53(. However, since the DAW 
frequency is much smaller that the dust charging frequency v%, one notices that Zdi and <j) 
are in phase. Subsequently, there appears a decrease of the DAW phase speed, as shown in 
the subsection below, 

From Eqs. (1), (2), (5) and (6) we readily obtain 



:|+^i)(v 2 -a£)-a> 2 



d 

4:7!-eZ d0 (— + ui)n d i, 



(9) 



where k D = (k 2 + kf) 1/2 = 1/X D , k e , t = l/\ De ,Di, andfc,= (47ra d n d0 ) 1/2 ■ 
Furthermore, from Eqs. (3) and (4) we have 



1 + T„ 



d_ 
dt 



.7^2 ~ V TdV > n dl + - 



Pd 



V 

dn dl 
dt 



(10) 



where Vrd = (/-^B^V^d) 1 ^ 2 is the effective dust thermal speed. Equations (9) and (10) are 
the desired equations governing the linear propagation of the modified DAWs. 

A. Modified DAWs 



Within the framework of a plane-wave approximation, assuming that n^i and </> are pro- 
portional to exp(— iut + zk ■ r), we Fourier analyze (9) and (10) and combine the resultant 
equations to obtain the dispersion relation for the modified DAWs 



l + k 2 X 2 D + "« 



K\ 2 D v 2 



(v\ — iuj) [uj{uj + iud) — k 2 V^ d + iujuj v } 



0. 



(11) 



where uj = k 2 Cj(l-R) > 0, C d = u pd X D , u v = r]*k 2 /(l-iur m ), and 77* = (£ + 4r]/3)/m d n d0 . 
Several comments are in order. First, for v d , <C \u\ <C 1/t to , fx, we have from (11) 



^ = ^ T^W^W D ' (12) 

where Uj = V^ d + rj^k 2 /r m . Second, in the long wavelength limit, viz. k 2 X 2 D <C 1, Eq. (12) 
reduces to 

^ = k ' U d + n , Tji , v (13) 
(1 + k 2 q \ 2 D v 2 jv x ) 

Third, the well-known frequency of the DAW, in the absence of dust grain correlations, dust 
fluid viscosities, dust-neutral collisions, polarization and ion pressure effects, and DCFs, can 



be obtained from Eq. (11). We have the famous result 43] 



kC d ft A\ 

- = (1 + ^)1/2- ( 14 ) 

Since C d = Z do [n do kBTi/n i0 m d (l + ai)} 1 ^ 2 , where at\ = n e0 Ti/nioT e , we observe from Eq. 
(14) that the phase speed (00 /k) of the long wavelength (in comparison with A_d) DAWs is 
proportional to [fcgTj/m^l + «i)] 1//2 , dictating that the restoring force in the DAW comes 
from the pressures of the inertialess Boltzmann distributed electron and ion fluids, while the 
dust mass provides the inertia to sustain the wave. The wave dispersion [the /c 2 A|,-term in 
Eq. (14)] arises owing to the departure from the quasi-neutrality condition in the perturbed 
density perturbations. 



B. Twisted DAWs 



We now discuss the possibility of a twisted DAV beam in the long-wavelength limit (viz. 
k 2 X 2 D <C 1), in which case the quasi-neutrality holds. Here we have 

n dl = -—^—(l + a 1 + a 2 )cj), (15) 

where «2 = cid^&B^i^do/ ^ie 2 n i0 . Equation (15) depicts that the dust density compression 
is possible, since 4> < (due to negative charge on dust grains) for the DAWs. 

In order to study the property of a TDAW, we consider the limit A|,V 2 <C <fi in Eq. (9) 
and combine the resultant equation with Eq. (10) to obtain the simple wave equation 

(V 2 + K 2 ) n dl = 0, (16) 
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where we have assumed that ridi is proportional to exp(— iut), and have denoted K 2 = u 2 /U 2 
with Ul = VS d + Cj(l - R)/(l + k 2 \ 2 D v 2 /vi). We have assumed that k 2 X 2 D < 1 and 
Vd,V*k 2 < \v\ < l/r m , 

We now seek a solution of Eq. (16) in the form 

n d i = ni(r)exp(ikz), (17) 

where n/(r) is a slowly varying function of z. Here r = (x 2 + y 2 ) 1 ^ 2 and is the propagation 
wave number along the axial (z—) direction. By using Eq. (17) one can write Eq. (16) in a 
paraxial approximation as 

2ik— + Vl )ni = 0, (18) 



dz 

where k = K = tu/U* and = (l/r)(d/dr)(rd/dr) + [l/r 2 )d 2 / d6 2 . We have used cylin- 
drical coordinates with r = (r, 9, z). 

The solution of Eq. (18) can be written as a superposition of Laguerre-Gaussian (LG) 
modes 60|, each of them representing a state of orbital angular momentum, characterized 
by the quantum number I, such that 



ni = n p iF p i(r, z) exp(U9), (19) 
P l 

where the mode structure function is 



F p i{r, z) = F pl X\ l \L^(X)exp(-X/2), (20) 

with X = r 2 /w 2 (z), and w(z) is the DA beam width. The normalization factor F p i and the 
associated Laguerre polynomial L p (x) are, respectively, 



1 P i 



and 



^l V2 , (2D 



L}[{x) = ^jP_ [xl+PeM _ x)l (22) 

where p and I are the radial and angular mode numbers of the DAW orbital angular mo- 
mentum state. In a special case with / = and p = 0, we have a Gaussian beam. 

The LG solutions, given by Eq. (19), describe the feature of a twisted DAV beam carrying 
OAM. In a twisted DAV beam, the wavefront would rotate around the beam's propagation 

8 



direction in a spiral that looks like fusilli pasta (or a bit like a DNA double helix), creating 
a vortex and leading to the DAV beam with zero intensity at its center. A twisted DAV 
beam can be created with the help of two oppositely propagating three-dimensional DAWs 
that are colliding in a dusty plasma. Twisting of the DAWs would occur because different 
sections of the wavefront would bounce off different steps, introducing a delay between the 
reflection of neighboring sections and, therefore, causing the wavefront to be twisted due to 
an entanglement of the wavefronts, and take on the shape of the reflector. Thus, due to 
angular symmetry, Noether theorem guards OAM conservation even for a longitudinal DAV 
beam. 

The energy flux of the DAV beam is given by W^Vn, where V s (= du/dk.) is the group 



velocity of the DAV beam and its energy density reads 



61] 



d IE 



2 



iy = -Mc,k)]_ fc ^L, (23) 

where E = ik(f>, and for u^, -C \lo\ <C l/V m ,Ti and k 2 \ 2 D C 1 we have e(a>, k) = {l/k 2 \ 2 D ) + 
(/CgZ/ 2 /fc 2 ^i) — ujp d (l — R)/{u) 2 — k 2 V^ d ). Equation (15) exhibits a relationship between the 
electrostatic potential^ and the dust density perturbation n d i- It turns out that the energy 
flux is independent of the mode number I. 

III. SUMMARY AND CONCLUSIONS 

To summarize, we have presented a new dispersion relation (11) for the DAWs in an 
unmagnetized dusty plasma that is composed of weakly correlated Boltzmann distributed 
electrons and ions in the wave potential, and strongly correlated highly- charged dust grains 
that follow the viscoelastic dust momentum equation. In our investigation, the effects of the 
polarization force and DCFs are incorporated. It is found that the contributions of both the 
polarization force and DCFs are to reduce the frequency of the DAWs. Furthermore, we have 
discussed the possibility of a twisted DAV beam. The latter can trap charged dust particles 
and transport them from one region to another. The present study of a twisted DAV beam 
can be useful for diagnostic purposes when the DAW frequencies are near the infra-sonic 
frequencies in laboratory, space and cosmic dusty plasmas. In closing, we mention that the 



importance o 



by Harwit 



OAM of electromagnetic waves in the astrophysical context was recognized 



62], and is also the subject of current interest 63j, |64( in connection with twisted 
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ultrasound pulses. 
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